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Abstract 

The forward-backward (FB) multiplicity correlations in two windows separated in 
'fr^"' rapidity and azimuth are analyzed in the framework of the model with independent 

^ ■ identical emitters (strings). Along with the short-range contribution, originating 

^ ! from the correlation between multiplicities produced by a single emitter, the long- 

range contribution, originating from the fluctuation in the number of strings, is taken 
I into account. The dependencies of the correlation coefficient on the rapidity and 

^ ■ azimuthal acceptance of the windows and on the corresponding gaps between them 

^ ■ are studied. It is shown that the analysis of these dependencies enables to separate 

, the contributions of two above mechanisms. We see that the presence of the pair 

^ i correlations between particles produced by one string turns the strings into non- 

^1 poissonian emitters. The connection of the FB multiplicity correlation coefficient 

with the double inclusive cross section, the two-particle correlation function and the 
r-| ■ so-called untriggered di-hadron correlation is also traced. The suitable observables 

are proposed. 

> . 

^ ■ 1 Introduction 

p : 

! For a long time, considerable attention is devoted to the experimental [l]-[6] and theoretical in- 
^ I vestigations of the forward-backward (FB) correlation - the correlation between multiplicities 
j—{ • Up and of charged particles produced in two separated in rapidity windows ("forward" 
^ ■ and "backward") in high-energy pp and A A collisions (see [7] -[11] and references therein). 
The main problem in the analysis of this correlation is the separation of the so-called "vol- 
ume" contribution, originating from the event-by-event ffuctuation in the number of emitting 
sources. 

In paper [17] it was suggested to use for this purpose the information on the event mul- 
tiplicity in an additional third rapidity window, but as discussed in [IB] it complicates the 
interpretation of obtained results. In present paper we argue that the investigation of the FB 
correlation between multiplicities in windows separated both in rapidity and in azimuth can 
enable not only to separate the volume contribution, originating from the ffuctuation in the 
number of sources, but also to obtain the important quantitative physical information on the 
magnitude of this ffuctuation in the processes under consideration. 

We also show that the traditional deffnition of the FB correlation coefficient leads to its 
strong dependence on the acceptance of the windows, with the correlation coefficient going 
to zero with the acceptance. As consequence the results obtained for the windows of different 
width can't be compared directly. In this connection we propose suitable observables for the 
FB correlation studies, which have some ffnite limit when the acceptance go to zero. 

To check our observations we use the simple two stage model [H [101 ED], inspired by a 
string picture of hadronic interactions. In this model one suggests that at the initial stage 
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of interaction some number N of strings are formed, which then are considered as identical 
independent emitters of observed charge particles. In previous note [21] we have considered 
only the long-range (LR) part of the correlation, originating from the fluctuation in the 
number of emitters (strings). 

In the present paper we also take into account the short-range (SR) correlation between 
particles produced by a single string. This SR correlation can arise due to very different 
physical processes such as the details of string break up, the formation and decay of clusters, 
resonances or minijets during the string fragmentation. We show that the presence of such 
SR correlation along with influence on the FB multiplicity correlation turns the string into 
non-poissonian emitter. 

The paper is organized as follows. In Sec. II we discuss the different version of the 
definition of the FB correlation coefficient and generalize these definition for the case of 
windows separated both in rapidity and azimuth. In Sec. Ill the connection of the FB 
correlation coefficient with two-particle correlation function is traced. On this base in Sec. IV 
we propose alternative suitable observables for the studies of the FB multiplicity correlations. 
In Sec. IV we shortly discuss the correspondence between the FB correlation and the so-called 
untriggered di-hadron correlation. 

In Sec. V we go to the model calculations and introduce the pair correlation function of 
a single string. In Sec. VI in the framework of the model we find the resulting expression 
for the FB correlation coefficient. In Sec. VII we analyze the typical dependence of the FB 
correlation coefficient on the distances between windows in rapidity and azimuth and discuss 
the separation of the contributions of two mentioned mechanisms. 

Appendix A describes the calculation of integrals over rapidity and azimuth windows. In 
Appendix B we present the alternative derivation of the basic formula for the FB correlation 
coefficient. 



2 Definition of the FB correlation coefficient 

Traditionally [H [21 HJ |5] the FB correlation coefficient is defined as a coefficient b in linear 
regression 

(^ij)n^ = a + 6n^ . (1) 

In this case 



(2) 



where D^^ is the variance of the multiplicity in the forward window 

Dn, = (nl) - {n^f . (3) 

Clear that the value of such defined correlation coefficient changes, if one will change 
independently the acceptances of the forward and/or backward windows. To avoid this trivial 
influence one can go from np and to the relative or scaled observables [22] i/p = Up/ (up) 
and z/^ = n^l (n^). In these observables {i^b)^^ — ^rei + ^rei 

In some papers [3l |6] the following symmetrized form of (|2]) is also used 

^ (npUs) - {np){nB) , . 
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for which one can prove that \bgy^\ < 1. Note that in the case of symmetric windows, when 
{up) = {n^) and = Dn^, all these definitions lead to the same result 



Kel = Kym = b ■ (6) 

In present paper we study the correlation between multiplicities Up and in windows 
separated both in rapidity and in azimuth. Denote by 6yp, S(pp and Sy^, 5ip^ the width of 
the forward and backward windows in rapidity and in azimuth, and by yp, (pp and y^, (p^ 
- the positions of the centers of the windows. We'll also use the following short notation for 
the acceptance of forward and backward windows 

5p = 6yp6(pp/2Tr , 5^ = 5?/^(5(y9^/27r . (7) 

By 

Vfb ^Vf-Vb ^ fPB = fF-fB (8) 

we denote the distance between the centers of the windows in rapidity and in azimuth. 

These variables are simply connected with the gaps y^^p and ip^^p between window in 
rapidity and in azimuth, we'll not use the last variables, but note that: 



Syp 5yg 6(pp 5ipg 

yFB = — + Vgap + — . ^FB = + ^gap + > (9) 

or for symmetric windows, when 6yp = 6y^ = 6y and 6ipp = 6ip^ = 6ip, 

Vfb = Vgap + . fFB = "Pgap + ^^P ■ (10) 

3 Connection with two-particle correlation function 

One can express the FB correlation coefficient through the two-particle correlation function 
^2(1/1, 1/2; — ^2)- For this, we have to introduce the pi(?/, ^p) and P2(l/i, V^i; 2/2, ^2) - the one- 
and two-particle densities of charge particles: 

pi{y, ^) = -i—r ' ^2(2/1, v'l; ?/2, ^^2) = — 1 — 1— • (n) 

dy dip dyi dipi dy2 d(p>2 

Then for the forward acceptance interval, at y G 5yp and ip G 5(p>p, we have [23] : 

dydippi{y,ip) = {np) , (12) 

/ dyidipi dy2d'p2P2{yi,'Pi;y2,y^2) = {np{np -1)) 

•f&Vp&Vp •J&yp&'Pp 

and the same for a backward window Sy^Sip^. Meanwhile at yi G 6yp, ipi G 6(p>p and y2 G 6y^, 
<P2 ^ S'Pb "we have 



dyidipi dy2d(p2p2iyuVi;y2,V2) = {npnB) ■ (13) 

Recall that (up) is an average multiplicity produced in the acceptance SypSipp. 
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The formulae (fT2|) and (fT3l) are the base for the experimental measurement of the one- 
and two-particle densities of charge particles pi{y,ip) and P2(?/i, V'l; 2/2, V^2)- For windows of 
small acceptance in rapidity and azimuth we have 

(rip) {''T'f^r) 

{np{np — 1)) 

P2(2/f> ^f; 2/f> V'f) = (^^y^^^^y ■ (15) 
Due to the rotation invariance in azimuth one has 

Pi(y>V5) = Pi(y)/2vr , p2{yu^uy2,^2) = p2{yuy2]Vi - ^2)/{'2Tcy . (16) 

Then we introduce two-particle correlation function 6*2(^1,^/2; V'l — V'2) by a standard way: 

^/ N p2iyi,y2;vi- ^2) ^ 

<^2 i/i,2/2;<^i -<^2 = f f 1- 17 

Pi{yi)Pi{y2) 

By ([I2])-([I7]) we have 

{npUs) - {nF){nB) = Ip^ , (18) 
Dn^ = {np) + Ipp, (19) 

where 

(^f) = ^ / ^2/ Pi (2/), (20) 
2vr Jsy^ 

Ip^ = {27iy^ dyidipi dy2d<^2Pi{yi)pi{y2)C2{yi,y2]^i- V2) , (21) 
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^FF = (27r) / dyidifii dy2dip2 pi{yi)pi{y2)C2{yi,y2; - ^2) ■ (22) 



What gives for the correlation coefficient 



^FB ^FB 



Dn„ (rip) + Ipp 



(23) 



Further simplification of the integrals is discussed in Appendix A. 

For example, in the case of FB windows separated only in rapidity (i.e. when S(pp 
Sip^ = 27r), by (11051) we have 



{rip) = / dyp.iy) , (24) 



hB= dyi dy2pi{yi)pi{y2)C2{yi,y2) , (25) 
Iff= I dyi c/2/2 Pi (2/i)pi (2/2)^2(2/1, 2/2) , (26) 

JSyp JSyp 

where 

1 r 

^2(2/1,2/2) = - / 6*2(2/1, 2/2; V^) (27) 



4 



and we used that 

^2(1/1, y2; -v) = C2(z/i, 1/2; ^) , C2iyi, 1/2; (p + ^nk) = (72(1/1, 1/2; v^) • (28) 

Note that if p2(l/i, V^i; Z/2, ¥'2) = Pi(z/i, V5i)Pi(l/2, ¥^2), then there is no correlation: C2 = 0, 
IpB = 0, Iff = and by Dn^ = M (see [23]). 

For windows, which are small both in rapidity and in azimuth, and within which one can 
consider 6*2(1/1, ?/2; fi — ^2) and Pi{y) to be constant, we have 

{^f) = pi{yF)h , {^b) = PiiyB)^B , (29) 

IfB = M{nB)C2{yF,yB','^FB) , (30) 

Iff = {npfC2{yF, yF', 0) , (31) 

Dn, = {nF)[l + {nF)C2iyF,yF-^0)] , (32) 

and 

h = ^^^^ h = ^^F)C2iyF,yB;fFB) /oox 

M i + {nF)C2{yF,yF;0) ' ^ ' 

Recall our short notations and ([H]). 

We see that the correlation coefficient (jl]), defined in scaled variables, still depends through 
(up) on the acceptance 6p of the forward window, that was observed earlier [20l [21] in a 
framework of a simple model. 

In the case when both small FB windows are situated in the central region, where one can 
suppose the translation invariance in rapidity: 

Piiy) = Po ' C2{yu 2/2; = 62(1/1 - 2/2; (34) 
the formulae (l29l) - (!33l) admit further simplification: 

(rip) = Pq5p , {rig) = p^S^ , (35) 

D„^ = K)[l + VoC2(0;0)] , (36) 

, _ , _ 5pPQC2{yFB\ ^Fb) /oyx 

Sb 1 + VoC2(0;0) • ^ ^ 



At last for large windows situated in the central rapidity region along with and ([35 
we must to use the formulae (fT9!) and (!23|l . with the following expressions for Ip^ and Ipp: 



^FB = Po(27r)M dyidifii dy2dip2C2{yi - y2\^i - ^2) , (38) 

Iff = Po(27r)"^ / dyidipi / dy2dip2 6*2(2/1 - 2/2; V'l - <y?2) ■ (39) 
The simplifications of the integrals in different cases are presented in Appendix A. 
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4 Alternative observables 

From (I29l)-(l33l) we see that for small symmetric windows at Sp = ^ we have 6 — )■ 0. 
This unpleasant dependence of the correlation coefficient on the width of the windows arises 
due to behavior of the variance in the denominator of (l23ll . Really by (!30|) and (132!) we 
see that in this limit Ip^ ~ SpS^ and ~ Sp. We can rid of this drawback if we normalize 
the correlator {ripn^) — {np){n^) by the product {np){n^) and introduce the observable 

{npUs) - {np){nB) Inp 



Then for windows, which are small both in rapidity and in azimuth, by fl30l) we have 

Kod = C2{yF.yB\'PFB) (41) 

or in the case of the FB windows, which are small only in rapidity and large i^^'^p = Sif^ = 2tt) 
in azimuth 

1 r 

TT Jo 

where we have take into account fl21|) - fl27|) . We see that in contrast with b the 6^^^ has a 
finite limit at small acceptances of windows. 

Another possibility, as it follows from ( !T9|) . is to use for the normalization instead of Dn^ 
and Dn^ the differences Dn^ — {rip) and Dn^ — {n^) and to introduce 

"rob — / I • l^-JJ 



Then again for windows, which are small both in rapidity and in azimuth, by (1301) and (l32l) 
we have 

\/C2{yp, yp] 0)C2{yB, ys'^ 0) 

or 

_ C2{ypB;VFB) 

~ C2{m ' ^ ^ 

if both of these windows are situated in the central rapidity region, where the translation 
invariance in rapidity takes place. 

In other case when the FB windows are 5ipp = 6(fB = 2vr in azimuth and small in rapidity, 
we have 

Kob = , ^^^'"'"^ , (46) 

where C2{yp,yB) is defined by f l27|) . In the central rapidity region C2{yp,yB) = C2{ypB) ^"^^ 
we have 

"--fir' (^^' 

We see that in contrast with h the h^^^ as the h^^^ has a finite limit at small acceptances of 
windows. 

Note that the definition fHSl) is closely connected with so-called robust variance [231 121 



Dn - (n) 

^ = I \2 ■ 4J 
6 



By (gOD and (03]) we have 

Kot = -7=== ■ (49) 

Emphasize that the traditionally defined ([2]) correlation coefficient b is also proportional 
to the two-particle correlation function C2{yp,yB', ^fb) (^^e (!33|) ). but the proportionality 
factor depends on the width of windows and goes to zero at 6p = 6^ 0. 



5 Untriggered di-hadron correlation 

The following alternative definition of two-particle correlation function is also in use 

C = S/B-1 , (50) 

where 

S - (51) 
dAy dAip 

Here Ay = yi — y2 and Aip = ipi — (p2 are the distances between two particles in rapidity and 
in azimuth, and one takes into account all possible pair combinations of particles produced 
in given event in some rapidity interval {Yi,Y2). The B is the same but for the case of 
uncorrelated particle production. 

At this definition in contrast with (fT7|) one implies from the very beginning that the 
translation invariance in rapidity takes place. Namely, that for any yi and y2 belonging to the 
interval {Yi, Y2) the result depends only on Ay = yi — y2- (AH the pairs with the same value 
of difference yi — 2/2 contribute to the same bin of the multiplicity distribution, irrespective of 
the value of {yi + 2/2)/2, see also the discussion in [IE].) This assumption is reasonable only 
in the central rapidity region at high energies. It means that we suppose that in the interval 
{Yi,Y2): 

piiy) = po ^ P2{yi^y2\v) = p2{yi-y2\v) (52) 

(see formula (!34|) ). 

In this case we have for the enumerator of fl50p : 

S{Ay,A(p)= dyidy2p2{yi-y2]Aip)5{yi-y2- Ay) (53) 

or in the case of commonly used symmetric interval {—Y/2,Y/2): 

S{Ay, A^) = p^{Ay; Aip) ty{Ay) (54) 

where the ty(A?/) is a "triangular" weight function fll02p . defined in Appendix A (see Fig 12]). 

In the denominator of fISUl) we should replace the P2{yi, 2/2; A^p) by the product Pi{yi)pi{y2), 
which due to the translation invariance in rapidity reduces simply to pg. Then 

B{Ay,A^)=pltY{Ay) . (55) 

Substituting into (15U1) we get 

C(Ay, A^) = P^^^y^^'P^ - 1 = C,{Ay, Av.) , (56) 
Po 
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where we have taken into account (1171) and (l34l) . We see that if the translation invariance in 
rapidity takes place within the interval (^1,^2), then the definition ( 150|) is equivalent to the 
standard one (IT7|) (see meanwhile the remark in the end of the section [7]). 

The drawback of this approach is that it supposes from the very beginning the translation 
invariance and hence can't be applied for an investigation of the multiplicity correlation at 
large rapidity distances, where the translation invariance in rapidity is not valid. At 
that by fl33l) . fj4T]) and (jH]) we see that the approaches based on the analysis of the standard 
([2]) or modified fHUl) . fH51) the FB correlation coefficients with two remote windows of small 
acceptance in rapidity and azimuth enable to measure the correlation strength C2{yi,y2] y^i — 
(P2) also in this case. 



6 The model. Pair correlation function of single string 

We now calculate the FB correlations in windows separated in rapidity and azimuth using 
the simple two stage model [H [TOl |20], inspired by a string picture of hadronic interactions. 
In this model we suggests that at the initial stage of interaction some number of strings 
are formed, which fiuctuates event-by-event with some variance -Dat = (A^^) — (A^)^ or scaled 
variance 

UN = Dn/{N) . (57) 

Note that the fiuctuation in the number of strings in pp and especially in AA collisions [28] 
is not poissonian and hence a; at 7^ 1. Its value depends on the collision energy. 

At next stage we consider these strings as identical independent emitters of observed 
charge particles. In previous note [2T] we have considered only the so-called long-range (LR) 
part of the correlation, originating from the fiuctuation in the number of strings. In the 
present paper we also take into account the short-range (SR) contribution, originating from 
the correlation between particles produced by a single string. 

To characterize the last property of the string we introduce, similarly to the consideration 
in the section [31 the two-particle correlation function for charged particles produced from a 
decay of a single string A{yi, 1/2; V'l — '^2)- For this purpose we at first introduce the Xi{y, (p) 
and \2{yi, V^i; ?/2, ^2) - the one- and two-particle densities of charge particles produced by one 
string. Then for given acceptance interval SypS(pp: 

dydipXi{y,ip) = {^p) , (58) 

I dyidipi c?Z/2'^V52 A2(yi,V5i;i/2,V52) = WW - 1)) 

and the same for backward window dy^dip^. Whereas 



dyidipi dy2d(p2X2iyu^i;y2,^2) = ifJ-Bf^p) ■ (59) 

SypSiPp JSygSipg 

The (lip) and (n^) are the average multiplicities produced by one string in the forward 6ypSipp 
and backward Sy^Sip^ windows. In general case due to the rotation invariance in azimuth 
one has 

= Ai(?/)/27r , A2(2/i,V2i;Z/2,</32) = A2(yi,y2;</?i -</?2)/(27r)2 . (60) 
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Similarly (fT7|) we can also introduce two-particle correlation function for charged particles 
produced from a decay of a single string by a standard way: 

A{yi,y2;'Pi-v>2) = — . . . . 1- 61 

Myi)My2) 



By dSHD-dHI]) we have 



where 



(/^b/^f) - = JpB > (62) 

= (^f) + JpF , (63) 

W) = ^ / dyXiiy) , (64) 



271 



5y 



F 



JpB = {2-n)'^j dyidipi ^1/2^992 Ai(yi)Ai(y2)A(yi, I/2; V^i - -^2) , (65) 

Jpp = {2ny^ dyidipi dy2dip2^i{.yi)^i{.y2)A{.yi,y2]'^i - '^2) ■ (66) 

JSypSipp JSypSipp 

Further simplification of the integrals are described in Appendix A. 

By fl63p we see that the presence of SR correlation turns the string into non-poissonian 
emitter. 

For small windows in which one can consider A{yi,y2] ^pi — 1^2) and Xi{y) to be constant 

(f^p) = Al(yF)^F . (/^b) = ^liyB)^B > (67) 

Jfb = h^B^i{yF)^i{yB)KyF^yB'^'PFB) > (68) 

<^ff = 4A?(2/f)A(i/f,2/f;0) , (69) 

where we have used our short notations 5p and 5^ ([8]) for the window acceptances. 

If the both small windows situated in the central rapidity region, where each string con- 
tributes to the particle production in the whole rapidity region, then due to the translation 
invariance in rapidity 

^i(i/) = /^o) Kyi^y2]^) = Kyi-y2]^) (70) 

and the formulae (I^7|) - fl69l) take the form 

(/^f) = /^O^F > (/^b) = l^O^B > (71) 

JpB = ^F^Bf^lMyFB'^ ^Fb) ' (72) 

Jff = 4/^oA(0; 0) . (73) 

Recall that yp^ and ^ps are the distances between the centers of forward and backward 
windows in rapidity and azimuth (|8]). 

Note that in the case of large windows situated in the central rapidity region along with 
( 17T|) one must use the formulae fl62|) and f l63|) with Jp^ and Jpp given by the following 
expressions 



Jfb = /^o(27r) ^ / dyidipi j dy2df2Hyi - 2/2; - ¥^2) , (74) 
Jff = /^o(27r)"^ / dyidipi / dy2d(p2A{yi - ^2; - ^2) (75) 



iSypSipp JSypSip 

(see Appendix A for further simplifications 
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7 The model. Resulting correlation strength. 



In a general case in the model with independent identical emitters [23]: 

pf(2/) = iVAi(y), (76) 

pUyi.y2\v) = iVA2(i/i,y2;</^) + iV(iV- l)Ai(i/i)Ai (1/2) • (77) 
Then the one- and two-particle densities of charge particles ffTTj) are given by 

pM = {p^y)) = {N)XM , (78) 

P2(yi,i/2;^) = {p^{yi,y2;v)) = {N)[X2{yi,y2;v) - Xi{yi)Xi{y2)] + (iV')Ai(i/i)Ai(i/2) , (79) 
and 

P2{yi,y2;^) - Pi{yi)pi{y2) = {N)[{X2{yi,y2;v) - Xi{yi)Xi{y2)] + D^Xi{yi)Xi{y2) > (80) 

where Dat is the event-by-event variance Dn = {N"^) — (N)'^ of the number of emitters. As a 

result we have the following expression of the two-particle correlation function 6*2(2/^, y^'i Vfb) 
(|T7|) through the pair correlation function of a single string A{yp,y^; (pp^) fl6T]) : 

^/ N + Hyi,y2;f) 

C2{yu 1/2; p) = , (81) 



where un is the event-by-event scaled variance un = Dn / {N) of the number of emitters ([5] 

In the case of FB windows which are small both in rapidity and in azimuth the flHTj) by 
(l33l) leads to the following formula for the FB correlation coefficient (jlj): 

h - ^^^K- ^^f)[^n + HyF,yB-^'^FB)]/{N) . . 

(n^) l + {np)[u:M + A{yp,yp;0)]/{N) ■ ^ ^ 

If these small FB windows are situated in the central rapidity region, where the translation 
invariance in rapidity takes place, then the A{yp,y^; Pp^) will depend only on the difference 
ypB = yp ~ yB of rapidities and the ( l82l) takes the following form 

6b l + 6pp,[uN + AiO;0)] ' ^ ^ 

where /Zq is the average rapidity density of the charged particles produced by one string. Note 
that in this case the basic formula (18T!) can be obtained also by an alternative way in the 
framework of the two stage model [9l [TOl [20] (see Appendix B). 

One can present the result for the FB correlation coefficient as the sum of two terms 
Kei = ^ref + ^^e?> where: 



l + (5^/ioK + A(0;0)] ' ^ ' 



and 



The first term depends on the acceptance 6p of the forward window, but doesn't depend 
on the distance between the centers of forward and backward windows in rapidity yp^ and 
in azimuth (fp^, which justifies the name of this contribution as the long range (LR) one. 
This contribution reveals itself as the common "plateau" when one plots the value of the 
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FB correlation coefficient 6 as a function of yp^ and ^pg- The level of this "plateau" is 
determined by the event-by-event fluctuation of the number of the strings (emitters) and 
can be used for the evaluation of the extent of this fluctuation. Note that at any fixed number 
of emitters there will be no such contribution, as cja? = D]y/{N) = 0. 

The second term is proportional to the pair correlation function A{yp^] (fp^) of a single 
string with some common factor depending on acceptance. In the plot of the FB correla- 
tion coefficient 6 as a function of yp^ and ipp^ this contribution manifests itself as some 
peaks above the level of the common "plateau" (see Fig |T] below). This justified the name of 
this contribution as the short range (SR) one. This different behavior enables separate the 
contribution of two mechanisms. 

We would like to emphasize that if the pair correlation function of a single string is equal 
to zero: A{yp^; ipp^) = 0, we still have nonzero FB correlation: 

"rel - ' ^^^^ 

1 + OpHqUJn 

which characterizes the event-by-event fluctuation of the number of strings A^. Note also that 
at A^yp^; (fp^) = by ( l63l) and ( l66l) the strings become a poissonian emitters and the answer 
(l86l) coincides with our result, obtained in [101 EOl El] for this case. 

For the alternative observables b^^^ ( HOjl and 6^^^ ( H3|l . introduced in section HJ in the 
framework of the model we have: 

UN + A{yp,yB;ippB) 

Omod - -jj^^ ' \P') 

UN + KVF^yB'^^FB) fc^a-. 
"rob = T-TT ^ • l^OJ 

UN + A{yp,yp;0) 

We see that they do not depend on the windows acceptances 6p, 6^ and have simple connection 
with the pair correlation function of a single string A^yp.y^; ^Pps)- 

In the case of large acceptance windows, within which one can't consider A[yi, y2] (p) to 
be constant, by ([I8])-(I23]) and (I76])-(IHI]) the formulae (IH2D, (EZD and (IHHD still can be used but 
with the following substitution: 

KyF^VB^^FB) ^ -, — n — r/ dyid^i 1 rfy2C^¥'2 Ai(?/i)Ai(y2)A(yi, 1/2; ¥'i-¥'2) , (89) 

where {^p) and (/i^) are the mean multiplicities, produced in the forward and backward 
windows by a single emitter (1M|) : 

= / dyXM , {^,p) = ^ / dyXM . (90) 

In the central region, due to translation invariance in rapidity, this substitution for large 
acceptance windows looks more simple: 



KVfb'^Vfb) ^ {.^yF^^F^yB^VB) ^ j dyidipl j dy2dip2 A{yi - y2]ipi - ip2) , (91) 
A(0; 0) {6yp6(pp)~^ / dyidipi / dy2dip2 A{yi - 7/2; ^1 - ^2) ■ (92) 

JSypSipp JSppSipp 
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Further simplifications of tlie integrals are discussed in Appendix A. 

For example, in the important case of symmetric {5yp = 6yp = 6y) FB windows, situated 
in the central rapidity region and separated only in rapidity (i.e. at 6ipp = Sip^ = 2tt 
acceptance in azimuth), by f llOip . f ll03p and fllOSp the formulae (PTj) and ( M2\i transform to 

Oy J-Sy 

A(0; 0) -> — / dyAiy)i5y - y) = 2 Aiy) dy - — ' yk{y) dy , (94) 



Jo Jo ^y 

where tsy{y) is a "triangular" weight function fll02p (see Figj2]) and 



A{y) = - [ d^A{y;^). (95) 

Jo 

We have used also that 

A{-y;ip) = A{y;^) , A{y; = A{y; ip) , A{y; ip + 2nk) = A{y; ip) . (96) 

In a conclusion of the section we note that if one uses the so-called di-hadron correlation 
approach, described above in section [5l for the experimental determination of the two-particle 
correlation function C{Ay,A(p) ( l50l) the result will depend on the details of "track and/or 
event mixing" used in that approach for the determination of B by the imitation of the 
"uncorrelated" particle production. 

If, as it was supposed in section [5] fl55|) : 

B{Ay,Aip) = / dyidy2pi{yi)pi{y2)S{yi-y2- Ay) = 

J-Y/2 

Y/2 

dy,dy2 (pf (z/i))(pf (1/2)) 5iy^ - y^ - Ay) = pi ty{Ay) = {Nfpl ty{Ay) , (97) 

-Y/2 

then by ([50]), (El) and (1791) with \i{y) = pg we get 

C(Ay, A^) = --^^(^V-^-P) . c,(Ay. A^) , (98) 

which under the assumption of the translation invariance in the central rapidity region cor- 
responds to the two-particle correlation function C2{yi,y2W)i defined by the standard way 
dUD, (compare with (IHT]) ). 

But if instead of (^71) one has 

/.y/2 

B{Ay, Aip) = / dy,dy2 (pf (l/i)pf (1/2)) S{y, - y2 - Ay) = {N')pl ty{Ay) , (99) 

J-Y/2 

as it frequently takes place in a di-hadron data analysis, then instead of (pH!) by (!50|) . (IMj) 
and (179|) we get 

C(Ay,A^) = ^A(Ay,Av.) , (100) 

which does not correspond to the standard two-particle correlation function 6*2(2/1, 1/2; V^), 
defined by f|T7j) . Compare (llOOp with fl98l) we see that in this case the resulting C{Ay, Aip) does 
not have an additional contribution reflecting the event-by-event fluctuation in the number of 
emitters. It depends only on the pair correlation function of a single string A(A?/, Ay?) and, 
therefore, is equal to zero in the absence of the pair correlation from one string. 
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Figure 1: The forward-backward (FB) correlation coefficient b (j5]) for symmetric windows of 
small acceptance 6yp = dy^ = 5y = 0.2 and 5ipp = 6ip^ = 6(p = 7r/4, calculated with taking 
into account the long-range (LR) and short-range (SR) contributions (1831) - (|85l) . at different 
values of a distance between the centers of the windows in azimuth (/? = (pp^ as a function of 
a distance between the centers of the windows in rapidity y = y^^ . 

8 Rapidity-azimuth dependence of the FB multiplicity 
correlation strength 

We analyze at ffist the FB multiplicity correlation strength in the most simple case of windows, 
which are small both in azimuth and in rapidity and situated in the central region, what is 
described by the formulae fl55]) - fl5S]) . 

As an illustration, the behavior of the correlation coefficient in the case of small FB 
windows, which are 6ipp = 6ip^ = 6ip = n/A in azimuth and 6yp = 6y^ = 6y = 0.2 in 
rapidity, with A{y; ip) expected from the Schwinger mechanism of a string break up are shown 
in Fig{T] [25] • (Note that in this case h^^i = b.) There are the large narrow nearside peak 

Vfb — Vfb — 0; some away- side structure at ipp^ = vr, which is smaller in 
amphtude and wider in rapidity, and the common "plateau" (pedestal) corresponding to the 
contribution of LR correlation, originating from the fluctuation in the number of emitting 
sources. 

Using the experimental data on the FB multiplicity correlation strength with small win- 
dows in azimuth and rapidity situated at different rapidity yp^ and azimuth ipp^ distances 
from each other, one can fix all parameters entering the SR (the parameters of the A{y; ip) func- 
tion) and the LR (the parameter un) contributions. (Note that only the products fioA{y; f) 
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and fio<^N are entering in the resulting formula (l83l) .) 

With the parameters, fixed in this way by the data on the FB correlation coefficient b 
with small acceptance windows, one can calculate (without any additional free parameters) 
using formulae fl9ip and (1921) the values of the FB correlation coefficient b for large acceptance 
windows, within which one cannot consider the A{y, ip) to be constant. 

For example, one can calculate the values of the FB correlation coefficient b in the prac- 
tically important case of symmetric windows separated only in rapidity, i.e. for 27r-windows 
in azimuth. In this case the resulting formula fl551) with the substitutions fl^5]) and fl^^ gives 
the dependencies of the FB correlation coefficient b on the width {5y) of windows and on the 
rapidity gap {ygap = Vfb ~ ^v) between them. 

9 Conclusions 

We have analyzed and compared the different deffnitions of the multiplicity correlation coeffi- 
cient. We see that the traditional definitions ([2]), (jlj) and ([5]) of the FB correlation coefficient 
lead to its strong dependence on the acceptance of the windows, with the correlation co- 
efficient going to zero with the acceptance. Hence, the results obtained for the windows of 
different width can't be compared directly. In this connection we propose suitable observables 
(HQ]) and (1431) for the FB correlation studies, which values have some nonzero limit when the 
acceptance go to zero. 

We have extended the definitions of the multiplicity correlation coefficient to the case of 
the windows separated both in rapidity and in azimuth. We have showed that this enables 
to separate the contribution, arising due to the event-by-event ffuctuation in the number of 
sources, and the contribution, originating from the pair correlation function of a single source. 
The last can reffect the very different physical processes such as the details of string break up, 
the formation and decay of clusters, resonances or minijets during the string fragmentation. 

In mid-rapidity region the first contribution, originating from the ffuctuation in the number 
of emitters, doesn't depend on the distance between the windows in rapidity and azimuth, 
which justifies the name of this contribution as the long range (LR) one. This LR contribution 
is proportional to the scaled event-by-event variance of the number of emitters un- 

The analysis of this contribution enables to obtain the important quantitative physical 
information on the magnitude of this ffuctuation in the processes under consideration. At 
that the strong non-linear dependence of the LR contribution on the window acceptances, 
taking place in the case of the traditionally deffned correlation coefficient ([2]) and (jl]), is fully 
specified (IHijl . 

The second contribution, originating from the correlation between multiplicities produced 
by a single emitter, depends on the distances between the centers of the backward and forward 
windows in rapidity and in azimuth ^Pp^. Its value ( 15^ is proportional to the pair 
correlation function of a single source A{yp^, (fp^) and decreases at large separations, which 
justifies the name of this contribution as the short range (SR) one. The dependence of the SR 
contribution on t/p^ and ipp^ can be retrieved from the experimental data on the multiplicity 
correlation with the FB windows, which are small both in rapidity and in azimuth. 

We also see that the presence of such SR correlation along with the inffuence on the FB 
multiplicity correlation inevitably turns a string into a non-poissonian emitter. 

We trace the connection of the FB correlation coefficient with the two-particle correlation 
function C2{yp,yB', Vfb) ^^nd show that the standardly defined two-particle correlation func- 
tion in general case is the sum of the contributions of the two above (LR and SR) mechanisms 
(see equation (IHTj) ). 



14 




y 



Figure 2: The "triangular" weight function arising due to phase space at integration over 
non-periodical FB windows (see Appendix A for details). 

We also discuss the so-called di-hadron correlation approach C{y,{p) fl50l) to the investi- 
gation of the multiplicity correlation function. The obvious disadvantage of this approach 
is that it supposes from the very beginning the translation invariance in rapidity and hence 
can't be applied for an investigation of the multiplicity correlation at large rapidity distances, 
where this invariance is not valid. At that we show that the approaches based on the analysis 
of the standard ([2]) or modified (140!) . (143!) the FB correlation coefficients with two remote win- 
dows of small acceptance in rapidity and azimuth enable to measure the correlation strength 

C2(2/_f5Z/b! Vfb) ^Iso in this case. 

Even in the mid-rapidity region, where the application of the di-hadron correlation ap- 
proach is justified, the results obtained by this approach depend on the details of "track 
and/or event mixing" used in this approach for the imitation of the " uncorrelated" particle 
production. As it was shown in the end of the section [TJ the obtained correlation coefficient 
C{y, if) can be equal to the standardly defined two-particle correlation function C-ziy, ^p) dHS]) 
and contain the sum of LR and SR contributions, or it can be proportional only to the SR 
contribution fllOOp in the dependence on details of the procedure applied. In the last case it 
looses the common "pedestal", the height of which is proportional to the scaled event-by-event 
variance of the number of emitters u^. 
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Appendix A. Calculation of integrals over rapidity and 
azimuth windows. 

For symmetric rapidity windows 6y^ = 6yp = 6y with the distance yp^ = yp — ys between 
their centers one has 




dyi dy2f{\yi-y2\)= dy f{\ypj^ + y\)tsy{y) , 



(101) 



where tsy{y) is a "triangular" weight function (see Fig|2]): 

tsyiy) = [9i-y)i5y + y) + e{y){5y - y)] e{5y - \y\) . 



(102) 
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The formula f llOip is valid for any distance between the centers of windows, in particular for 
coinciding windows. In the last case = and we have 



r r rSy rSy 

/ dyr dy2f{\yi-y2\)= dy f {\y\) Uy{y) = 2 dy f {\y\){5y - y) (103) 

JSyp JSyp J-Sy JO 

The same formula 

r r^v 
d^i j d^2f{\Vi-^2\)= j d^f{\^pB + ^\)ts^{^) (104) 

SiPp JStfig J —Sip 

is valid for the integration over azimuthal windows, but in this case one has also to take into 
account the periodicity: /(Iv^l) = f{\^ + 2Tik\). The last leads to significant simplification of 
the formula (11041) in the case of windows of full 27r acceptance in azimuth: 

27r /'27r nix 

dipf{\ippB + ip\)t2A^) = 2n dipf{\ip\) = An dipf{\ip\). (105) 

-2n Jo Jo 

So for large symmetric windows in the central rapidity region by fll0ip - fll03l) the formulae 
(|9T|) and fl92|) in general case can be written in the following form 

/Sy i>5ip 
dy dipA{ypB + y;^FB + v)tsy{y)tspi^), (106) 
-Sy J —Sip 

rSy rSip 

A(0; 0) ^ A{5y5^)-^ / dy / d^ A{y- ^) {5y - y) {5^ - ^) . (107) 
Jo Jo 

The 6y and 6ip are the width of the observation windows in rapidity and in azimuth, and the 
ypp and ipp^ are the corresponding distances between their centers. We imply that A{y; (p) 
satisfies the conditions ( l96l) . The same simplifications take place for the integrals Ip^ ( l38l) 
and Ipp (139!) in section [31 

Appendix B. Connection of resulting correlator and vari- 
ance with the ones for single emitter 

In the framework of the two stage model [H [TOl [20] in the case of the FB windows in the 
central rapidity region one can express the observable correlator (npn^) — {np){n^) and the 
variance Dnp through the correlator {fJ^pfi^) ~ (/^f) if^s) variance D^^ for one emitter: 

(npUj^) - {np){n^) = (A^)((/i^/i^) - {fip){fiB)) + Dn{iip){iib) , (108) 

Dn^ = {N)D,^+DN{fip)\ (109) 

where Dn = {N'^) — (N)'^ and (N) are the event-by-event variance and the mean number of 
emitters (see [21], [29] for derivation). These formulae are obtained for the windows situated 
in the central rapidity region, where each string contributes to the particle production in the 
whole rapidity range and hence the translation invariance in rapidity takes place. In this 
region by ([351) and ([71]) we have also 

(up) = 6pPo , (fip) = 6pfiQ , {hb) = SbPo , ifis) = ^Bf^o > Po = (^)/^o • (HO) 
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Recall that /Xq is the average rapidity density of the charged particles produced by one string. 

In the case of FB windows which are small in rapidity and azimuth by (fT8|) and (130|) we 
have 

^2[yF^yB^VFB) = 7 — w — V • (HI) 

{Up) (n^) 

Similarly by ( 162|) and ( 172|) we have the same for the pair correlation function Kixj^^yp] fps) 
of a single string: 

MVb, Vf: ^fb) = ^^^^y^;;;^^ ■ (112) 

Then combining the formulae f ll08p -( fTT2l) we again get the formula ( IHTl) of the text in this 
particular case. 
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